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A  Statistical  Derivation  of  the  Average  Degree  of  Polymerization 

in  a  Stirred  Tank  Reactor 

Neil  Dotson 

Department  of  Chemical  Engineering  and  Materials  Science 
University  of  Minnesota 
Minneapolis,  MN  55455 


Abstract 


TTie  statistical  approach  in  its  many  forms  has  been  a  powerful  tool  in  describing  the 
structural  evolution  of  nonlinear  polymers,  but  has  been  limited  in  its  application  to  batch 
processes.  This  restriction  is  removed  by  considering  the  paradigm  of  Af  homopolymeriration 
in  a  homogeneous  stirred  tank  reactor.  The  number-  and  weight-average  degrees  of 


polymerization  are  derived  from  recursive  arguments  which  properly  account  for  the  non- 
random  combination  of  monomers  caused  by  the  residence  time  distribution.  The  results 
in  agreement  with  those  derived  from  the  kinetic  equations.  This  route  is  not  necessarily  a 
simpler  or  more  powerful  approach  than  solution  of  the  kinetic  equations,  but  it  does  provide 
new  insights  into  this  particular  problem,  as  well  as  demonstrating  the  versatility  of  the 
statistical  approach. 
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Introduction 


Stepwise  polymers  produced  in  a  homogeneous  continuous  stirred  tank  reactor* 
(HCSTR)  can  exhibit  significantly  broader  distributions  than  their  counterparts  produced  in 
batch  reactors.  The  homopolymerization  of  A2  monomers  serves  to  illustrate  this.  In  a  batch 
reactor  the  molecular  weight  distribution  is  the  most  probable  (or  geometric)  distribution:* 

Pn  =  pn-*(l-p)  (1) 


where  Pj^  is  the  molar  fraction  of  n-mers,  and  p  is  the  extent  of  reaction  (the  probability  that  a 
functional  group  has  reacted).  This  distribution  gives  the  following  number-  and  weight- 
average  degrees  of  polymerization: 


DP„  = 


1  -  p 


DP 

1  -  p 


(2) 


(3) 


The  polydispersity  of  the  polymer  thus  has  an  upper  bound  of  two.  The  conversion  is  related 
to  time  t  through  the  Damkohler  number  Da  by  the  following  relationship: 


Da 

1  +  Da 


(4) 


where  Da  =  ktfcjo,  k  is  the  rate  constant  for  functional  group  (ii£il  monomer)  reaction,  f  is  the 
functionality  of  the  monomer  (here  2),  and  Cjq  is  the  initial  monomer  concentration. 

In  an  HCSTR,  on  the  other  hand,  the  distribution  takes  a  more  complicated  form:* 
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P  (2n  -  2)!  p"-^ 

^n"n!(n  -  l)!(i  +  p)2n-l 


(5) 


which  yields  the  following  average  degrees  of  polymerization: 


DP„  = 


n“  1  -  p 


1  +  p2 


DPu.= 


(6) 

(7) 


Here,  in  contrast  to  the  previous  case,  the  polydispersity  has  no  upper  bound.  The  conversion 
is  related  to  the  relative  rates  of  exit  and  reaction  through  the  Damkohler  number  by  the 
following  relation: 


P  = 


Vl  +  4Da  -  1 
I5i 


(8) 


where  Da  =  kSfcjo,  Cjg  is  the  entering  monomer  concentration,  and  the  mean  residence  time  0 
is  the  ratio  of  the  volume  of  the  reactor  to  the  volumetric  throughput  rate. 

Since  for  linear  polymerizations  broad  distributions  arc  easily  obtained  in  an  HCSTR, 
for  nonlinear  polymerizations  one  might  expect  the  situation  to  be  even  more  severe,  such  that 
the  gel  point  could  be  advanced  conn  ’e.  ably.  Several  researchers  have  investigated  the 
behavior  of  nonlinear  polymerizations  in  HCSTRs.  With  regard  to  step  polymerizations,  the 
first  work  is  that  of  Cozewith  et  al.,7  which  solved  both  the  steady  state  and  start-up  problems 
for  homopolymeiization.  For  finite  functionality  systems,  Gupta  and  coworkers^'^  have 
analyzed  HCSTRs  in  series,  the  effect  of  oscillations  in  feed,  and  the  effects  of  intramolecular 
reaction.  Hendriks  and  Ziff^  have  analyzed  the  problem  starting  with  the  Smoluchowski 
coagulation  equation,  and  have  looked  at  different  coagulation  kernels,  analyzing  both  A2  and 


4 


A«  homopolymerization  (which  describes  as  well  the  critical  behavior  of  the  Af  system). 
Nonlinear  free-radical  systems  have  also  been  studied.”^*^^ 

Common  to  all  of  these  analyses  is  the  use  of  the  kinetic  equations  describing 
polymerization.  This  approach  would  seem  necessary  because  the  various  statistical  methods, 
which  are  so  useful  for  describing  nonlinear  polymerizations,  have  been  held  (usually 
implicitly)  to  be  valid  only  for  batch  polymerizations.  ^  2,1 3  The  non-geometric  character  of  the 
distribution  of  eq  5  would  thus  seem  to  be  beyond  statistical  reach,  but  yet  Biesenberger  and 
Tadmor  justify  it  with  combinatoric  reasoning.  Thus,  this  alleged  restriction  is  not 
absolute.  Nonetheless,  the  question  must  be  asked;  how  are  the  ideality  assumptions  (equal 
reactivity  and  independent  reaction)  violated,  so  that  simple  application  of  statistical  methods  is 
invalid?  Then,  is  there  a  correct  statistical  -  ;  to  the  average  properties?  The  paradigm  of 

random  step  homopolymerization  of  f-functional  monomers,  analyzed  by  the  recursive 
method,^^  serves  to  answer  both  questions,  with  the  answer  to  the  latter  being  "yes." 


Statistical  (Recursive)  Derivation. 


DP,. 

Consider  the  homopolymerization  of  f-functional  monomers  in  an  HCSTR  operating  at 
steady  state,  for  which  cyclization  is  prohibited,  and  for  which  only  monomer  enters  the 
reactor.  To  derive  the  weight-average  degree  of  polymerization,  we  first  grab  a  (zeroth 
generation)  monomer  at  random,  so  that  we  can  write: 

DP,=  WfE(NA®"*)  (9) 


where  E(Na°“*)  is  the  expected  number  of  attached  monomers  looking  out  from  a  functional 
group.*5  We  must  realize  two  things:  (1)  the  probability  of  any  of  the  f  functional  groups  on 
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the  monomer  having  reacted  is  dependent  upon  the  residence  time  of  that  monomer,  and  (2)  the 
only  correlation  between  the  states  of  the  different  functional  groujjs  is  the  fact  that  they  all 
share  the  same  residence  time.  These  two  points  correspond  to  violations  of  the  equal 
reactivity  assumption  and  the  independent  reactivity  assumption,  respectively.  The  source  of 
these  nonidealities,  however,  is  not  intrinsic  to  the  reaction  chemistry,  but  rather  to  the  reactor 
configuration. 

As  an  HCSTR  is  perfectly  mixed,  the  residence  time  distribution  is  described  by  an 
exponential  distribution: 


p(t)dt-i  e-^®dt 


no) 


Given  a  residence  time  t  of  the  zeroth  generation  monomer  (see  Figure  1),  the  probability  of  a 
functional  group  having  reacted  is  given  by  the  following: 


p(t)  =  1  -  e-^^ 


(11) 


where  t  is  the  lime  scale  for  reaction,  equal  to  [k(l-p)fciol'l.  such  that  6/x  =  Da  (1-p).  The 
overall  conversion  is  thus  given  by: 


P  =  Jp(')P(.)d>  (12) 

which  is  consistent  with  eq  8. 

Conditioning  upon  residence  time,  we  can  then  write: 


E(NAO“0  =  jE(NAO“‘|t)p(t)  dt 


(13) 
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To  evaluate  E(Ny^o**t  1 1 ),  we  must  first  condition  upon  the  time  t'  at  which  the  reaction 
occurred  (see  Figure  1): 

t 

1 1 )  =  1  -  e‘^^  +  (f-1)  Je(N;^o“*  I  one  group  reacted  at  time  t')  -e‘(*'^  dt* 

(14) 

where  the  weighting  function  in  the  integran'^  ■  the  probability  density  of  reaction  at  time  t' 
given  a  residence  time  t,  which  is  consistent  with  eq  1 1.  To  close  the  recursion,  we  must 
further  condition  upon  the  residence  time  t”  of  the  first  generation  monomer  (see  Figure  1); 

E(Na°“‘  i  t )  =  1  -  +  (f-1)  j  fE(NAO“‘  1 1”)  p*(t"-t’)  dt"e‘’^^  dt'  (15) 

The  recursion  thus  does  not  close  algebraically,  as  it  does  in  the  case  of  a  batch  reactor,^^  but 
rather  gives  a  double  integral  equation. 

The  distribution  p'(t"-t’),  which  is  the  probability  density  that  a  monomer  has  a 
residence  time  t"  given  that  one  of  its  groups  reacted  at  t',  is  a  function  only  of  (l"-t’)  because 
the  residence  time  distribution  of  a  monc  v. ,  one  group  of  which  reacted  at  time  t',  is 
insensitive  to  the  fact  that  it  must  have  remained  in  the  reactor  for  a  subsequent  period  t'.  The 
distribution  is  given  by  the  following: 

p'(t)  dt  =  (—-<-  —  )  e'^®  dt  (16) 

0  X 

i.e.  simply  the  (properly  normalized)  product  of  the  probabilities  of  a  functional  group  escaping 
both  reaction  and  exit 

The  double  integral  equation  can  now  be  transformed  to  a  second-order  ordinary 


differential  equation,  and  the  result  is: 
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E"(Nj,ovt  1 1)  .  ^E'(N^out  1 1)  +  (  i+  1)  (  ECN^out  1 1)  +  ^  )  =  Q  (17) 

The  character  of  the  solution  will  differ  depending  on  whether  f=2  or  f>2;  this  is  to  be 
expected,  since  a  linear  system  is  very  different  from  a  nonlinear  one.  Before  proceeding  to 
the  solutions,  though,  we  can  answer  the  first  question  posed  in  the  Introduction. 

If  we  had  naively  assumed,  in  eq  l4,  that  the  first  generation  monomer  were  also 
randomly  chosen,  we  would  have  written  that 

1  one  group  reacted  at  time  t')  =  E(Ny^o“0  (18) 

The  recursion  would  thus  have  closed  algebraically,  and  we  would  have  obtained: 

oo 

1  +  fp(t)  p(t)  dt 

DPw  = - ^ -  (19) 

1  -  (f-1)  Jp(t)  p(t)  dt 


which  by  eq  12  is  simply  the  familiar  batch  reactor  resultl^  Thus  we  see  that  it  is  not  merely 
the  residence  time  distribution  that  is  the  source  of  the  nonideality,  but  rather  that  the  residence 
time  distribution  forces  non-random  combination  of  monomers.  That  is,  a  monomer  with  a 
residence  time  t'  does  not  combine  with  other  monomers  with  a  probability  proportional  to  p(t); 
rather,  the  residence  time  distribution  of  the  first  generation  is  different  from  the  overall 
distribution.  In  fact,  the  expected  residence  time  of  a  first  generation  monomer  is  0(1- 
(x/(0-rx))2),  which  is  less  than  the  value  for  the  zeroth  generation  monomer,  6.  The 
"skewing"  of  the  residence  time  distribution  will  also  occur  for  every  successive  generation, 
and  so  the  nonideality  may  be  called  a  long-range  effect  The  problem  can  nonetheless  be 
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solved  exactly,  as  shown  below.  (All  results  will  be  expressed  in  terms  of  the  conversion  p; 
see  eq  8). 

I.  The  Case  f  =  2. 

Under  the  initial  condition  ECN^o***  1 0)=0,  and  the  restriction  that  must  be  finite, 
eq  17  is  easily  solved  to  yield: 

E(NAO“Mt)=  -(-+  ~)t  (20) 

-t  t  6 


The  expected  weight  on  a  monomer  thus  grows  linearly  with  residence  time.  Integrating 
according  to  eq  13,  we  find  the  weight-average  degree  of  polymerization  to  be: 


as  before. 


DPw 


_  1 

"(1  -  p)2 


(7) 


n.  The  Case  f  >  2. 

The  solution  of  eq  17  depends  upon  whether  we  are  operating  below,  above,  or  at  the 
critical  conversion  Pc,  or  the  cotresponding  critical  Damkohler  number,  Da^.: 


2(f-2)(-y/  i  *rh  • ' ) 

(21) 

°^"4(f  -  2) 

(22) 

as  found  from  the  discriminant  of  the  characteristic  equation.  We  will  examine  each  case  in 
turn. 

A.  p<pc.  With  the  initial  condition  that  E(Ny^o“*  1 0)=0,  and  the  condition  that  as  p  (or 
Da)  approaches  zero  DP^  must  approach  unity,  the  solution  is: 
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E(SAO^i\t)  =  j^(t^^-\)  (23) 

where 

^  =  )  (24) 

20  \  T  T 

In  contrast  to  the  linear  A2  homopolymerization,  for  the  nonlinear  case  the  expected  weight  on 
a  monomer  grows  exponentially  with  residence  time.  The  weight-average  degree  of 
polymerization  is  given  by: 

DPw=(jr:^^((l-p)W(l-p)2-4(f-2)p  )  -j4t  P5) 

in  agreement  with  the  solution  which  can  be  obtained  from  the  kinetic  equations  (see 
Appendix), 

B.  p=pc.  Here  we  obtain: 


(26) 


which  in  turn  gives  DP^r  as  follows: 

DPw  =  2f^f4  (^“^^ 

Thus  the  largest  DP^^  to  be  obtained  is  four,  for  f=3;  as  f  approaches  infinity,  the  largest  DP^^, 
goes  to  two.2  Thus,  despite  the  broad  distributions  ,  we  do  not  obtain  high  polymer.^  (The 
polydispersity  here  is  deceptively  low,  for  at  any  Damkohler  number  all  moments  above  a 
critical  moment  are  divergent.^) 


C.  p>P(;.  '^or  this  case,  the  roots  of  the  characteristic  equation  are  complex  and  thus 


the  solutio’-  is  of  the  form: 


E(NaO“‘ It)  =  e‘^2e{cjsin(^'>^ 4(f-2)  ®  (1  +  Vl  ) 


(28) 


C2Cos{^^j^f-2)UlA■l  )  }  -  fTT 


an  equation  which  makes  no  physical  sense,  as  it  allows  E(N^oui  |  be  negative.  Indeed, 
DPy,  is  found  to  be  given  by  eq  25,  and  thus  itself  to  have  an  imaginary  component.  We  thus 
conclude  either  that  steady  state  operation  is  not  possible  for  p>Pc,2  or  that  the  steady  state  is  a 
gelling  one.6  Practically,  of  course,  gelation  would  make  operation  of  an  HCSTR  impossible 
anyway,  such  that  from  an  engineering  viewpoint  p^  gives  the  critical  conversion  (and  hence 
Da^.)  above  which  steady  state  operation  is  not  possible. 


DP„. 

The  number-average  degree  of  polymerization  can  always  be  calculated  from  simple 
stoichiometric  reasoning,  being  related  only  to  the  number  of  reacted  groups  and  not  to  the 
details  of  connectivity.  Thus: 


"  1  -  (f/2)p 


(29) 


in  agreement  with  the  result  from  the  kinetic  equations  (see  Appendix). 


Discussion 


For  the  sake  of  illustration,  we  examine  the  case  f=3;  all  of  the  conclusions  drawn, 
however,  are  independent  of  f.  The  steady  state  values  of  the  number-  and  weight-average 
degrees  of  polymerization  are  shown  in  Figure  2,  which  emphasizes  that  high  polymer  is 
simply  not  obtained  in  a  single  HCSTR.  However,  while  is  finite  (and  in  this  case  equal 
to  four)  at  the  critical  conversion,  it  should  be  noted  that  (dDF^^/dp)  is  infinite.  Insight  can 
also  be  gained  from  comparison  of  the  HCSTR  results  with  those  from  a  batch  reactor;!^  we 
c  in  compare  on  the  basis  of  either  conversion  or  Damkohler  number.  For  DP^,  Figures  3  and 
4  show  these  comparisons;  on  either  basis,  the  HCSTR  product  exhibits  a  higher  DP^  than  the 
batch  reactor  product,  as  expected.  This  is  obviously  not  due  to  increased  conversion  at  a 
given  Da;  indeed.  Figure  5  shows  that  the  DPn  for  the  HCSTR  product  is  lower  than  that  from 
a  batch  reactor  at  the  same  Da,  indicating  a  correspondingly  lower  conversion.  (A  comparison 
of  DPn  on  the  basis  of  conversion  is  not  shown,  because  there  is  no  difference  between  the  two 
on  this  basis.) 

In  conclusion,  the  recursive  method  has  been  extended  to  HCSTRs,  proper  care  being 
taken  to  account  for  residence  time  distribution  and  the  accompanying  non-random  reaction. 
The  results  obtained  are  in  agreement  with  kinetic  derivations,  and  although  this  derivation  is 
not  necessarily  simpler  or  more  powerful  than  the  kinetic  derivation  (indeed,  we  have  only 
solved  the  steady-state  situation),  it  does  offer  insight  into  how  the  structure  changes  with 
residence  time  of  the  monomer.  More  importantly,  it  demonstrates  the  power  of  the  statistical 
method  to  account,  rigorously,  for  nonideality,  in  this  case  a  violation  of  both  the  independent 
and  equal  reactivity  assumptions  enforced  by  the  reactor  configuration.  While  the  ideas 
presented  here  could  be  extended  to  other  more  practical  chemical  systems  (e.g.  Af  +  B2)  or  to 
other  types  of  reactors,  perhaps  they  will  find  greater  utility  in  the  analysis  of  nonidealities 
intrinsic  to  the  reaction  chemistry. 
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Appendix:  Derivation  of  DP^  and  DP^  from  kinetic  equations  for  f>2. 

The  kinetic  equation  describing  the  polymerization  is  essentially  Smoluchowski's 
coagulation  equation  with  source  and  sink  terms,®  and  with  a  kernel  appropriate  to  Af 
homopoly  meri  zation :  ^  ® 


^  +  =1  X  [(^2)i+2][(f-2)j+2]CiCj.k[(f-2)n+2]Cn£K^^^  (Al) 

®  i+j=n  i=l 


where  Cno  is  the  entering  n-mer  concentration.  The  first  three  moments,  (i=0,l  ,2),  of  the 

distribution  evolve  according  to  the  following  expressions;!? 


^  .  k  {  (l/2)(f-2)Vi^  +  2(f-2)mMo  +  2Mo^  } 


dt  0 


^  k  {(f-2)2^^'  '»)P2lll  +  4pi2} 


V.  here  pjo  entering  i*  moment.  The  concentration  of  unreacted  functional  groups  [A] 
(which  equals  {(f-2)|ii+2po})  obeys  the  following  equation: 


For  steady-state  operation  (i.e.  the  time  derivatives  vanish)  and  only  monomer  entering 
(Mo(^l^lO=^^2(^‘^lo)»  solution  of  the  above  (algebraic)  equations  gives: 


Mo  = 


if,  Vl+4Da  -  1 
*  2  ^  ■  2Da 


)) 


^10 


(A6) 


=  cio 


(A7) 


^2  = 


/'fVl-4(f-2)Da  _2_^ 

2Da(f-2)2  ■  f-2  / 


(A8) 


From  eq  A5,  the  conversion  is  found  to  be  given  by  eq  8,  and  upon  defining  DP^=P2/^i  and 
DPn=pi/}iQ,  eqs  25  and  29  are  obtained. 
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